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Abstract 

The leading divergences of the generating functional for Green functions of quark currents 
between one-nucleon states are calculated with heat kernel techniques. The results allow for a 
chiral invariant renormalization of all two-nucleon Green functions of the pion-nucleon system 
to 0{p^) in the low-energy expansion. 
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1. The modern treatment of the pion-nucleon interaction as an effective field theory of the 
standard model was pioneered by Gasser, Sainio and Svarc who applied the methods of chiral 
perturbation theory (CHPT) to the vrA^ system. However, the chiral expansion with baryons 
is hampered by the presence of the nucleon mass, which stays finite in the chiral limit. It was 
then shown by Jenkins and Manohar that the methods of heavy quark effective theory |Q 
allow for a systematic low-energy expansion of baryonic Green functions in complete analogy 
to the meson sector. 

The purpose of this letter is to perform the complete renormalization of Green functions 
of quark currents between one~nucleon states to 0{p^) in the chiral expansion. This requires 
a chiral invariant calculation of the divergent part of the corresponding generating functional 
at the one-loop level. Many specific Green functions of the vrA^ system have already been 
analyzed to one-loop accuracy ^, ^, ^. However, only the full divergence structure of 0{p^) 
presented here permits a complete renormalization of all Green functions: nucleon form factors, 
vrA^ ^ vr . . . irN, 7*iV ^ vr . . . vrA^, W*N n . . . irN, etc. Results will be given for two hght 
flavours only. All details and the extension to chiral SU{3) will be presented elsewhere. 

We start from QCD with two light flavours u, d coupled to external hermitian fields 0: 

C = rl^cD + 9-7^ (V,. + + 75^m) 9 - q{S - il^P)q^ " ( d ) • 

5, P are general 2-dimensional matrix fields, the isotriplet vector and axial- vector fields V^,A^ 
are traceless and the isosinglet vector field is included to generate the electromagnetic current. 
At the effective level of pions and nucleons, V* couples directly only to nucleons since the pions 
have zero baryon number. 

Explicit chiral symmetry breaking is implemented by setting S = ^A = diag {rriu^md)- The 
chiral group G = SU{2)l x SU{2)ji breaks spontaneously to SU{2)v (isospin). It is realized 
non-linearly Q on the Goldstone pion fields 0: 

a{4>) ^ gLa{<t>)h{g,<pr\ g = igL,9R)eG 

Cr{<P) ^ gRiR{4>)h{g,(l))'\ (2) 

where Cl,(,r are elements of the chiral coset space SU{2)l x SU{2)ii/ SU{2)v and the com- 
pensator field h{g, cj)) is in SU{2)v- The more familiar quantity f7(</>) = C-R('/')Ci('/')^ transforms 
linearly under G. In the standard "gauge" with = =: n we have U = v? . 
The nucleon doublet ^ transforms as 

* = ( ^ ) -^' = %,'^)* (3) 

under chiral transformations. The local nature of this transformation requires a connection 

T/. = - ir,)U + did, - ie,)^} (4) 

in the presence of external gauge fields 



V = + A^,, = - A, (5) 
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to define a covariant derivative 

v,^ = id, + r^-ivp^. (6) 

The effective chiral Lagrangian for Green functions with at most two nucleons is |jl], |^ 



C2 = —{D^UD''U^ + xU^ + X^U) = —{u^u'' + x+) (8) 



^eS = C,M + ^MB (7) 

Cm = C2 + J~-i + . . . 

X = 2Bo{S + iP), x± = u^XU^ ^ux^u 

l-MB — i-^jv + -^-ttAT + ^vrAT + • ■ ■ 

611 = (9) 

where m,gA are the nucleon mass and the neutron decay constant in the chiral limit and (. . .) 
stands for the trace in flavour space. 

The generating functional of Green functions Z[j,r],ff\ is defined by the path integral 

^iZ[j,v,n] = [dud^d^ exp[i{5M + Smb + j d^x{fj^ + *r/)}]. (10) 

The action Sm + S]\j b corresponds to the effective Lagrangian (0), the external fields V, A, S, P 
are denoted collectively as j and ij, f/ are fermionic sources. 

2. Heavy baryon CHPT 0] can be viewed as a clever choice of variables for performing the 
fermionic path integral in (|lO|). By shifting the dependence on the nucleon mass m from the 
nucleon propagator to the vertices of the effective Lagrangian, the integration over the new 
fermionic variables produces a systematic low-energy expansion. 
In terms of the velocity-dependent fields N^,Hy defined as 

Nv{x) = exp[imv ■ x]P^"^{x) (11) 
Hv{x) = explimv ■ x\P~ {x) 

the pion-nucleon action Smb takes the form 

Smb = j d^x{N,AN, + H„BN, + N.^'^B^^'^H, - H,CH,} (12) 
A = ivV + qaS ■u + A^2) + A(^-i) + ... 

B = i J/^ - iAy . + 5(2) + 5(3) + . . . 

C = 2m + ivV + qaS ■ u + C(2) + C(3) + . . . 
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In A and C, the only dependence on Dirac matrices is through the spin matrix 

S''='-^5a^^''v,, S-v = 0, ^' = -^1- (13) 

Rewriting also the source term in (|^) in terms of N^,H^ with corresponding sources 

= e^'^^-^P+r/, = e^'^^-^p-r/ , (14) 

one can now integrate out the "heavy" components to obtain a non-local action in the 
"light" fields |9|, |l^, 11]. At this point, the crucial approximation of heavy baryon CHPT 
is made: the action is written as a series of local actions with increasing chiral dimensions by 
expanding in a power series in 1/m: 



With this approximation, the integration over A^^, reduces again to completing a square (the 
fermion determinant is trivial to any finite order in 1/m) with the final result 

gj^Ii,'?.*?] — ]V J ^fly^j^i{SM+ZMB[u,j,Pv,Rv]) ^IQ-'j 

where 

+ R^C-^B{A + 7°pt^°C-^S)-S°^H°<^"^^.; - RvC-^Rv}- (17) 

The functional integral ([lO| ) has been reduced to the mesonic integral (1T6|). From here on, the 
standard procedure of CHPT can be applied: the action in the functional integral (16) is 
expanded around the classical solution Uc\[j\ of the lowest-order equation of motion. Integration 
over the quantum fluctuations gives rise to a well-behaved low-energy expansion for Z[j, r], fj] 
like in the meson sector. 

For u = Uci, we can read off the tree-level nucleon propagator from the right-hand side of 
Eq. ( p!7[ ) (see Fig. 1 of Ref. Q for an artistic impression of this object): 

^mI [j,V,v] = Zmb [uc\ [j\ ,j,Pv,Rv]- (18) 

By construction, Z^fj^ is independent of the arbitrary vector v. It is instructive to check this 
independence for the free propagator corresponding to Q 

A = ivd, B = i^^, C = ivd + 2m. (19) 

Inserting these operators into the right-hand side of Eq. (17) and reexpressing p^, R^ in terms 
of the original fermionic sources, one finds indeed the free massive fermion propagator 

Zji7|[0,7?,7?] = / d'^xfi{x){i ^ + m){a + m^)-'^r]{x). (20) 



On the other hand, any given order in the chiral expansion of ZMBijjV^v] '^ill general not 
be independent of v because a change in v involves different chiral orders (reparametrization 
invariance [^]). This is the price one has to pay for a systematic low-energy expansion. 

^C~^ is not expanded in this case. 
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3. We now turn to the calculation of ZmbUj^^v] up to and including 0{p^). Since we only 
consider nucleons (rather than antinucleons), we can drop the sources to the order considered. 
To 0{p^), ZjvfB[j, r/, r/] is a pure tree-level functional. The relevant pion-nucleon Lagrangian is 

A(^i) =iv - V + qaS ■ u 5(1) =iJ7^- ■ 



A(2) = Ci{u-u) + C2{{v-uy) + C3X+ + C^{x+) (22) 

where -Fl,h are the field strengths associated with the isotriplet gauge fields £,r in (^). Since 
the loop contributions set in at O(p^) only, neither gA nor the coupling constants Ci, . . . , C7 of 
0{p^) will have to be renormalized. 

To calculate the loop functional of O(p^), we expand 

C2 + Ca- pt,A^^)p^ (23) 



in the functional integral (16) around the classical solution Ucib]- A convenient choice of fluc- 
tuation variables ^ is given by ||^ 

iRi^) = n(0ci)e*«('^^/', = nt(<Ad)e-^«('^)/', = t iO = 0, e(0ci) = 0. 

(24) 

Following Gasser, Sainio and Svarc we expand £2 to 0{^^), £4 to 0(0 and A^^^ to 0(^2) 
to arrive at the diagrams of Figs. 1,2. The sum of the reducible diagrams in Fig. 2 is finite and 
scale independent 0. For the irreducible diagrams of Fig. 1, we need 

^(1) = iv -V + qaS ■ u 

= iv ■ Vci + qaS • + - [u • Uci, (,] - gAS ■ VdC 

+ '-^v- Vci e + yK, [S • nei,^]] + Oif) (25) 

From now on, the index cl will be dropped. All mesonic quantities are to be taken at the classical 
solution of the lowest-order equation of motion. 

The diagrams of Fig. 1 correspond to the following generating functional: 



Zhv[j,Pv] = J <fx(fx'(fy(fy'pv{x)A^^^{x,y)[T.i{y,y')5^{y - y') + S2(y,y')] • ^ij:){y\x)pv{x) 

(26) 
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where A^^ is the propagator for Ny in the presence of external fields. The self-energy functionals 
Si, S2 are given by 

^i{y,y') = -^{in[Gij{y,y')v- d'jk -V ■ dijGjk{y,y')]Tk + gA[Ti,[S ■ u,Tj]]Gij{y,y')}{27) 
^2{y,y') = -^Vi{y)Gij{y,y')Al^^{y,y')V,{y') (28) 

The differential operator d^- in 3-dimensional tangent space is related to the previously defined 
covariant derivative as 

V^^ = ^r,</„ e = ^r,e., d^ = S^,^'^ + l^, ll", = -\{^''[n,r,]) . (29) 



It acts on the meson propagator Gij: 

G = {d^d^^ + a)~\ aij = l{[Uf„Ti][Tj,u'']+x+{Ti,Tj}). (30) 

o 

4. The self-energy functionals T,i{y,y'), T,2{y,y') are divergent for y' y. In order to extract 
the divergences in a chiral invariant manner, we use the heat kernel representation of propagators 



(see Ref. [13] for a review) in d-dimensional Euclidean space. The divergences will appear as 
simple poles in e = |(4 — d). The corresponding residua are local polynomials in the fields of 
0{p'^) and can easily be transformed back to Minkowski space. 

In Euclidean space with d dimensions, the inverse of the elliptic second-order differential 
operator 

D2 = -d^d^ + a (31) 

can be constructed as an integral 

poo 

G(x, y) = D^\x, y) = / d\G{x, y; A) (32) 
over the heat kernel G{x,y; A) with its asymptotic expansion for A ^ 0+ : 



Gix,y;X) = j^e-^^-yr'/^'Y.-nix,y)X- 

^ + D2^G{x,y;X)=0, G{x,y;0) = d'^ix - y). (33) 

The differential equation for G{x, y; A) yields recursion relations for the Seeley-DeWitt coeffi- 
cients an, which are matrices in tangent space of 0(p^"). Since the divergences of Si,S2 are 
0{p^), we only need the following derivatives of oq, ai in the coincidence limit x ^ y, which is 
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to be understood in the list below: 
ao = 1 



df^dyGQ = d^ao d„ = a„ = lim [d^anix, y) - a„(x, yh^{y)] 



dxdf,di,ao = 2dxdf,ao d^, = ^{dx^if^y + d^^^ixu) I^lv = d^^i, - d^j^ + bfiHu] (34) 

ai = -a dx'jfj.u = dxi^iu + b\,l^lu] 
df,ai = ^d^'jfj^i, - ^dfj_a 
ai d^ = -^du^i^u - ^d^a. 
In d dimensions, we have m, m 

ai{x,x) _4-d 

]imd'^^{dy)G{x,y) = ^—^md;{dy)ai{x,y). (35) 



Inserting the coincidence relations for ai into Si in (pTj), performing the summation over the 
tangent space indices k and transforming back to Minkowski space, one obtains the divergent 
part of Si, which is manifestly of 0{p^): 

sr(y,y) = j^,^^(y) (36) 



^i(y) = --^^'"'^f^uv" + — {{S ■ u,u ■ u + x+} + S ■ u{u ■ u + x+ 
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+ 2u^{u^'S -u) - {S -uu-u) - {S ■ UX+)) (37) 



5. The divergences of T,2{y, y') are due to the singular behaviour for y ^ y' oi the product of 
(derivatives of) propagators 

G,,{y,y')A-^^^{y,y') . (38) 

In accordance with locality, the divergences are of the generic form (with an overall chiral 
dimension 3) 

^^^^F(y, y')D(y; v)5\y - y') (39) 

with a field monomial F{y, y') and a differential operator D{y; v) of at most third order. 

The missing ingredient is the heat kernel representation for the nucleon propagator A.'^iy 
Below, only results needed for the present analysis are collected. However, the following repre- 
sentation will clearly also be useful in heavy quark effective theory. 
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Since ^(i) is not an elliptic differential operator, we write 

A^-^^=iV\VV^y^, V = iA(^i) (40) 

and set up a heat kernel representation for the inverse of the positive definite hermitian operator 
A = VV^, where 

T> = —v-V — igAS-u + rjv'^, r/ — > 0+ 

A = VV^ = -(v - Vf + R + T]^ (41) 
R = igA{v -yS ■ u) - g\{S ■ uf. 

To facilitate the transformation back to Minkowski space, it is preferable to keep track of the 
dependence on v'^ in Euclidean space. To get a well-defined representation for A~^ without 
infrared divergences, one must keep r/ 7^ in intermediate steps. The heat kernel representation 
for A~^ can be motivated by remembering that A is essentially a one-dimensional differential 
operator in the direction of v. 

A-\x,y) = / dTE{x,y;T) (42) 
Jo 

-^ + A^E{x,y-r)=0, E{x,y;Q) = 5\x - y) 



E{x,y;T) = — exp 
\/47rT 



v-{x- y)]2 2 
r] T 



4r 



g{x -y;v)^bn{x,y)T'' 



9{x; v)= j j^^^ik ■ v)e-"'-^ vdg = 0. 

The differential equation for E produces again recursion relations for the coefficients bn, which 
are matrices in 2-dimensional flavour space. The following values of the bn{x,x) are needed for 
the present analysis: 

60 = 1, (^;.V)"6o = 
bi = -R, V ■Vbi = -^v - VR. (43) 

In addition to the Seeley-DeWitt coefficients a„, 6„ and their derivatives, we need the coin- 
cidence limits of (derivatives of) products of the meson and nucleon propagator functions 



„2 , 



Gn{x)= d\ Jt^ , Em[x;v) = g{x;v) dr . (44) 

Jo (47rA)^/^ Jo V4vrr 

The corresponding singularities can best be extracted in Euclidean cZ-dimensional Fourier space 



14|. The following list contains all singular products appearing in S2: 



Go{x)Eo{x-v) I 5\x) e = \{A-d) (45) 

Goix)vdEo{x;v) ~ ^ „ v Od'^ix) (46) 



S^S,d^d,Go{x)Eo{x;v) ~ .^^S^S,6^,{vdf6''{x) (47) 

2 



S^S,d^d,Go{x)vdEo{x;v) ~ ~ 3(4^)2/ ^'^^ " ^)''^''(^) (48) 
S^S,d^d,Go{x)Ei{x;v) ~ V^''!^) (49) 



5M'5.5^^-Go(x)t;-a^i(x;f) ~ .^^Sf,S,6^,v ■ d6\x) (50) 

(47rj^e 

5;.5.9;.9.Gi(x)^o(^;^^) ~ -r^r^S^S,6^,6''{x) (51) 



5^5,a^a,Gi(x)t;-9ii;o(x;7;) ~ -j-4^S^S,6^,vdS''{x). (52) 

The actual calculation of the divergent part of S2 is rather tedious. The final result can be 
expressed in terms of a local functional S2(y) defined in analogy with (^) : 

^r{y,y') = (4;^^2(y)5''(y - y')- (53) 

Transforming back to Minkowski space, one finds 

5^2(2/) = i ' ~'~ ((^ ' '")^)]'^ • V + -V • m(v • Vv ■ u) + -(w • u{v ■ Vv ■ u)) 

+ QA |~^^ ■ ^("S" -uv -u) + ^(S • ti(u • uf) - S^v''\r^y,v ■ 



3, . „ 3, „ — 1 

32 



-(u^(i;Vn^)) - -(w • V-u • n) + - — • 9xh 



- ^[(^^•^i) + 2n-n + ^(x+)]?;-v| 

+ si {-^5 • 7x(z; • V)2 - 25^^(5 • nr^,)5^ _ i(^ . V5 • n)^; • V 
_ . v)25 . ^.) - ^u^iu^S .u) + ^{S-u{u-u + x+)) 

- ^ {X+,S ■u} + ■u{x+)^ 

+ I [2(5 • uf - 4((S • uf)]v ■V + ^{v-VS-u)S-u 

+ ^5 • u{v -VS-u)- 4(5 • u{v ■ VS ■ n))| 5^ 

+ 9AS,[l{S-uf-^{{S-uf)^S'^. (54) 

For the purpose of renormalization, this lengthy expression for S2 is not in the most suitable 
form . For the following counterterm Lagrangian, powers of the spin matrix are reduced with 
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the relations 

{S^, S"} = ^{vf^v" - g^"), [S^, 5^ = ie'^^'P^'vpS^. (55) 
In addition to various SU{2) identities, I have also used the curvature relation 

1 i 

^fiu = ^[u^^u^,]- -f+^^. (56) 

6. With the conventions of Ref. [Q] for separating the finite part in dimensional regularization, 
we decompose Zirrij^Pv] in ( p6| ) into a finite and a divergent part, both depending on the 
arbitrary scale parameter /i: 



^i{y,y')5\y-y') + ^2iy,y') = sf(y, y'; ;u)54(y - y') + (y, y; m) 

U\y-y')[^,{y) + ^2{y)] (57) 



2A(/i) ,4, 



F2 

,d-4 



The generating functional Z[j,p^] can now be renormalized by introducing the counterterm 
Lagrangian 

4?(^) = J^:^J2B^N,{x)0,{x)N,ix) (58) 

with dimensionless coupling constants Bi and field monomials Oj(x) of 0{p'^). In analogy to 
Eq. (|57|), the low-energy constants Bi are decomposed as 

A = i?[(M) + (47r)2/3iA(/i). (59) 

The /3j are dimensionless functions of qa designed to cancel the divergences of the one-loop 
functional. In Table 1, the operators Oi are listed together with the coefficients /3j. 

This completes the renormalization of Green functions to 0{p^). The sum of the irreducible 
one-loop functional Zij-j-[j, py] and of the counterterm functional is finite and scale independent: 

Zi„[j,Pv] + (4^^7)2 Yl^i J d'^^d'^^'d^yMx)A{i]ix,y)Oi{y)A~Yy,x')pt,{x') = 
d^xd'x'd'ydS'p,{x)Al^^{x, y)[Sf(y, y'; fi)6\y - y') + T.f{y, y'; p) 
+ j^5\y-y')Y.Bl{pmy)]A-^{y\x')p.{x'). (60) 

Since the same is true for the sum of the reducible contributions of Fig. 2, the total generating 
functional of 0{p^) has been rendered finite and scale independent. 

The renormalized low-energy constants BI (p) are measurable quantities satisfying the renor- 
malization group equations 

P-^BUp) = -A (61) 

implying 

B:ip) = B:{po)-(3.log^. (62) 

Mo 
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Table 1: Counterterms and their /^-functions as defined in Eqs. |59| ) 



i 


Oi 


A 


1 




-(l + 55i)/12 

(4 - 5a)/8 
gAi^-g\)/8 
gA{6-6gl + g\)/l2 
-9a{8 - g\)/8 


2 


i[u^,V^'v ■ u] 


3 


i[v ■ u,v • Vv • u] 


4 


S ■ u{u ■ u) 


5 


Uf^{u^S ■ u) 


6 


S-u{{vuf) 


7 


V ■ u{S ■ uv ■ u) 


-9^/12 
(l + 55i)/24 
5A(4-5i)/8 
-(l + 55i)/6 


8 


[X-,v- u] 


9 


S ■ u{x+) 


10 




11 


iS^^v'^lf+^.^v ■ u] 


9A 


12 




-g\{4 + 3g\)/16 
-9\l^ 


13 


vxe^^''PSp{{v ■ Vup)u,) 


14 




-9\l^ 
-^9\ 


15 


i{v ■ 


16 


V- V S ■ uv ■ V 


9\ 


17 


{u ■ u)iv ■ V+ h.c. 


-3(7^(4 + 3(7i)/16 
(8 + 95i)/16 


18 


{{v ■ u)^)iv ■ V+ h.c. 


19 


(u • V5 • u)v ■ V+ h.c. 


-95a/16 
-5i(4 + 5i)/4 
9\ 


20 


• V+ h.c. 


21 


vxe^i^^PSpUpU^v • V+ h.c. 


22 


ivxe^'^'PSpUp^v ■ V+ h.c. 
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The operators in Table 1 constitute a minimal set for renormalizing the irreducible self- 
energy functional for off-shell nucleons. As long as one is interested in Green functions with 
on-shell nucleons only, the number of operators in ( |58[ ) can be further reduced by invoking the 
equations of motion for the nucleons. We shall not pursue this option here as it requires some 
additional discussion. 

Partial results for one-loop divergences in chiral SU{2) have been obtained before. In par- 
ticular, the results in Table 1 agree with those of Ref. ||l^. Some of the coefficients have also 
been calculated recently by Park, Min and Rho: the coefficients q in Eq. (43) of their paper Q 
agree with Table 1, except for C4 which should be a factor two smaller. 

7. In addition to phenomenological applications for chiral SU{2), the methods of this paper 
can be applied and extended in various ways: 

i) Chiral S[/(3). 

ii) Inclusion of higher baryon states, especially in view of recent results on the 1/Nc expansion 
for baryons (see Ref. |l^ and references therein). 

iii) Renormalization of the one-loop functional of 0{p'^). 

iv) Non-leptonic weak interactions of baryons. 

v) Combination of heavy quark and chiral symmetries. 
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Figure Captions 

Fig. 1: Irreducible one-loop diagrams. The full (dashed) lines denote the nucleon (meson) 
propagators. The double lines indicate that the propagators (as well as the vertices) have 
the full tree-level structure attached to them as functionals of the external fields. 



Fig. 2: Reducible diagrams of O(p^). The cross denotes counterterms from the meson La- 
grangian £4 of O(p^). 
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